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1. SUMMARY 


This report is a statement of progress on NASA Grant NSG-3048 during 
the six month period from September 1, 1979 to February 29» 1980. During 
this period > the researches at the University of Notre Dame were directed 
by Professor Michael K. Sain; and the funded research assistant was Mr. 
Stephen Yurkovlch. Mr. R. Michael Schafer continued studies developing 
out of researches under this grant; his support, however, was drawn from 
fellowship funds made available by the University of Notre Dame. Over- 
seeing technical aspects of the grant at Lewis Research Center was Dr. 
Kurt Seldner. 

The major emphasis of this status report lies In the continuation of 
nonlinear modeling researches Involving the use of tensor analysis. Pro- 
gress has been achieved by extending the studies to the controlled equa- 
tion 

X - f(x,u) 

and by considering more complex situations. Included herein are calcula- 
tions Illustrating the modeling methodology for cases In which x and u 
take values In real spaces of dimension up to three, and In which the de- 
gree of tensor term retention Is as high as three. The quality of the 
controlled tensor models has been most encouraging; and preparations are 
now under way to begin applications to the QCSEE digital simulation. 

Though not funded by Grant NSG-3048 during this period, certain lin- 
ear multivariable studies growing out of earlier grant work are described 
briefly, for completeness. 
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2. STATUS OF CURRENT RESEARCH 


This section reviews the status of grant researches carried out during 
the slx>month period from September 1, 1979 to February 29 > 1980. 


In the previous Semi-Annual Status Report, for the period from March 
1, 1979 to August 31, 1979, the initial studies on nonlinear modeling by 
means of tensor ideas were presented. Basically, those results were con- 
cerned with homogeneous systems of the form 

X • f(x), 

for X a member of some finite dimensional real vector space X. If 

f(0) 5* 0, 

then the equilibrium point x^ satisfying 

f(x^) - 0 

can be translated to the origin. Without loss then. It Is assumed that 

f(0) - 0. 

Then f(x) is to be understood in terms of its power series expansion 
about the origin, with each real component 

f^ : X -► R 

of the function 


f : X -*• X 

leading to series terms of the form 


X. X. ••• X. . 

*^1 2 ‘'p 

Such terms, though nonlinear, are nonetheless p-llnear functions 

R** -*• R, 

which by basic tensor theory can be expressed in terms of linear functions 
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on the tonsor product (p times) of R with Itself. In the aggregate, the 
collection of such linear functions can be understood In terms of a matrix 
operator on the tensor product (p times) of X with Itself. This matrix 
operator can be derived. If f la known, or determined empirically if f 
Is not. The preceding Seml^Aimual Status Report considered both viewpoints. 


For control problems, of course. Interest would center on an equation 
of the form 


X - f(x,u), 

where u Is a member of another finite dimensional real vector space U. 
The pair 

(x,u) 

is then an element of the product space 

X X U, 

which is sometimes denoted 


X0U. 

The focus of this report Is on this case, and the details are presented In 
Section 2.1 following. 
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All discussion and examples to this point have concerned systems with 
homogeneous nonlinear vector differential equation representations. While 
emphasis has been centered around the linearization of nonlinear systems 
based on identification of the total operator L, the real thrust of the 
concept has been on identification of the individual linear operators as- 
sociated vlth each censor product term. It is this vein that Che discus- 
sion will now pursue, with the inclusion of control (or input) variables in 
Che system. Certainly Che individual linear operators can be obtained from 
Che partitions of Che large matrix operator if so desired. Let it suffice | 

to say chat either could be considered a "by-product" of the other. 

The initial venture of this chapter is Co establish Che notation and 
basic structure of the necessary tensor algebra for inclusion of forcing 
functions in the system. With this framework secured, the existing tech- 
nique is adapted Co accommodate nonhomogeneous multivariable systems; here 
a concept of basis ordering is introduced. To exhibit the application of 
Che ideas developed, several examples are treated and resulting simulations 
plotted. First, two second-order systems are examined— that is, with two 
states and two controls. Progressing from these, a three state, three 
control example is Inspected. Throughout, the choices for the forcing fune 
cions are sine wave inputs, with amplitude and frequency varied for each 
input. Again, intuition falls in these nonlinear equations as Co Che 


choice of Initial conditions; therefore, saell initlel conditions 

ere chosen to ensure complience with the feasible operating region of each 

system. 

To distinguish between the concepts in this chapter and those developed 
previously, the usual symbols for the state vector and control vector will 
be employed: x and u, respectively. 


2.1.1 Formulation: Forcing Functions Included 

The aim of this section is to extend the concepts treated thus far 
(the fflultlllnearlty of the higher degree nonlinear terns, for example) to 
functions of two variables, namely states and controls. With the notation 
convention established here, the software adaptation and representative ex- 
amples can be discussed in the sequel. 


Consider the general set of nonlinear vector differential equations 

X • f(x,u) 


where the state vector x is an element of R and the input vector u 


is an element of R 




the Hessian matrices 


Define the Jacobian matrices 
lax^ij ’ 

ij ^ ’ 


3^f. 


3x 


p l*lp2p««e ptl 


as 


2 2 
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and Ch« croaa darlvatlvt matrices 


r 

i: 

[ 


l!!i 

dxdu 


* i " l»2,«»«,n 


by 


2 2 


With a fixed operating point (x^^Uq), the general equation may be written 

as 


f(x,u) ■ f(XQ,U^) 


^T- 

ix 


if 


3u 






(x-x^) 


(u-Uq) 


n 

+ 1/2 I e,(x-xj’ 1 

i-1 ^ ^ 3x 

n 3^f 

* 1/2 Z 

i«l 3u 


(Xq.Uq) 


(x-Xq) 


n 3^f 


(xQtUo^ 


(u-u^) 




(u-Uq) 


+ . . 


where (') denotes transposition and e^ is the real vector of zeros 
except for the i-th entry which is unity. 


Suppose now tnai. 


ir~ '<v"o>- 


so that by the substitutions 
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ic follows that 


• If “ * ■ • 

Thsrtfors, ths original diffarantlal aquation la now eonvartad to tha fol- 
lowing slmpllflad fora: 

. df| 


3x 




9u 


(Xq.Uq) 


n a^f 

+ 1/2 I a. z’ 1 

1-1 ^ 3x 

n aV 

+ 1/2 I a. w' y 

1-1 ^ au^ 


n a^f 

+ I a' 

1-1 ^ 






(Xq.Uq) 


+ . . . 


To Introduce tha notions of the tensor algebra, consider for Instance 
the term 


n a^f 

I *1 

1-1 ^ ax 




and regard it as a bilinear function 

n a^f 

I *1 O’ -4 

1-1 ^ ax 


( ). 




Because of the bilinearity, this function can be expressed in terms of 


th« tensor product ss s linear function operating on s • z. Call this 


particular linear operator L 2 Q so that 

2 
n 


1/2 I e z' 

1-1 ^ 3x* 


(Xo.Uq) 


z - 1/2 L 2 Q (z ® z). 


The notation for the linear operator Is such that the first subscript cor- 
responds to the first vector, z, and the second subscript corresponds to 


w. Thus, by this convention, 

lii 

3x3u 


I ‘i 


1-1 


w - 1^2 (z • w) . 


(’‘O'V 


NOW define 


ho- If 


L 

01 3u 




(Xo.Uq) 


as the usual "A" and "B" matrices of a standard linearization 


z - Az + Bw. 


With this established, the differential equation can now be expressed as 

* • 4 o * ^ 1 '' 

+ 1/2 Ljo (* • t) 1/2 Lo2 (w • tf) 

+ (z ® w) . . . 

The general term In such an expansion can be represented as 

Ljj^ (z ® z ® . . . ® z ® w ® w ® ® w) , 

J times k times 
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to thoe, floally, 

* • I I TTrr (8 •...•* ® w ® ...« w) . 

J «0 k -0 ■> > 

j elatt k tiatt 

Agalof th« aotivation for tueh a darivacion it tha face that all of 
tha art lloaar oparatort. Tha nuabar of thaat oparatora for a glvan 

ascpantlon dapanda, of eouraa, on cha daalrad accuracy of tha approxlaacion. 
Contidar» for axaapla, cha tabla pratanCad in Figurt 2.1 for cha cypical 
aytcaa ccnaitcing of cha tcaca vaccor x and Cha concrol vaccor u. 


linaar oparatort 


4 o **01 
4o 4l ^02 
ho hi hi hi 
ho hi ^22 hi ^04 
^50 hi hi ^21 ^14 ^05 


camt rtCainad in approxiaacion 
tcandard linaarixacion 

aacond dagraa 

chird dagraa 
fourch dagraa 
fifch dagraa 


Figure 2.1 

Accordingly; for an approxlaacion rtcaining up Co and Including all 
Chird dagraa Cantor carat , nine Individual naad ba IdanClflad; for 

fourch dagraa carat, cha nuabar trould be fourctan, and to on. 

In cha raaalning tacclont of Chit chapCar cha IdanClficacion achaaa 
it datcrlbad and llluacracad. Howevar, a alnor polnC thould ba aanClonad 
hara bafora procaadlng fur char. Vhan cha ara IdanClfiad, cha pra- 
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I 

I 

1 

j 

1 

1 

« 

r 

.» 

I 


ccdint constant tarn aasoeiatad with aach>»that ia» 1/jlkl— la aupprasaad 

Into tha operator itself. This obviously will cause no loss of 

since any one of the constants aay be extracted (if so desired) from the 

corresponding operator by a simple matrix sealing operation. Furthermore, 
* 

as will be seen in the next section, the operators that are identified 
are of reduced size as compared to those used in the discussion to this 
point. 


2.1.2 Software Adaptation 

The major difference between the algorithm here and the one employed 
in the first and second order linearizations of the previous two chapters 
is, cf course, the necessary addition of control variables. However, a 
change more fundamental to the basic technique is the need now for a 
greater number of time points in the sampling process. This is basically 
due to the behavior of the input forcing function, taken to be sinusoid 
for the examples treated in following sections. Obviously, the sampling 
period must be somewhat smaller than the period of oscillation of the in- 
put. For ease in programming, a constant sampling rate is used. Another 
major adaptation in the Software is the development of a scheme for logical 
ordering of the terms which arise from the tensor products. 


Consider again a series expansion, given now in terms of the vectors 
X and u, 

* ■ ‘■10 * * ‘•01 “ 

* ‘■20 * * ‘•11 ‘•02 “ 

4’ X • X • X 4- . . . 
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i 






Since the tesk is to identify the parameters of the only, the new vector 


equation 


X - [Lj_q-: Lqj_; L2q1 Lq2: I-jq: • • • 1 


X 

u 

X • X 
X 9 u 
u 9 u 
X ft X 9 xl 


would be constructed. 

Recall that for the case of the homogeneous system several redundant 
terms appeared in the vector stacked with tensor product terms. In such a 
case it was argued that those terms, with the corresponding columns of the 
identified operator, could be eliminated for purposes of the identification. 
Such is the case here; for Instance, consider the product term 

X ft X ft u 

which corresponds to the L 2 j^ operator. For three states and two Inputs, 
this product can be constructed according to the previous scheme: 

xftxftu^xft (xftu) 

f«. 

X ® 




X ® 


x.u. 

X 

1 

c 

1 1 

i ; 

X»u. 

x^u, 

2 1 

2 : 


x^u 

3 1 

3 i 


Again, the term (x ® u) can be considered as a slx-dlnenslonal object, and 
stacked In one (row) vector. Thus, 

Tx, 


X ® X O u 


'"l"! 


*l“2 


*2^1 


*2^2 


XjUi 


*jUjl 


r 2 

*i“i 

2 

*l“2 

*l*2“l 

X1X2U2 

X1X3U1 

X1X3U2 


*2*l“l 

X2X1U2 

2 

X2U1 

2 

*2“2 

*2*3“l 

X2X3U2 

f 

X3X1U1 

Vl“2 

X3X2U1 

X3X2U2 

2 

X3U1 

2 

X3U2 



^n 18-dlmenslonal object. Therefore, the linear operator associated 

with this product term would have dimension 3 x 18. But notice that there 
are six Identical terms in the product: 

*l*2"l ■ ’‘2’'l'*l’ 

’‘l*2“2 * ’‘2*l'*2’ 

XiXjUi - XjXj^Uj^, 

X1X3U2 - XjXj^U^, 

X2X3U1 - *3’'2“l* 

X2X3U2 - X3X2U2. 

Eliminating the redundant terms of x 0 x ® u would leave a vector of length 
12. Likewise, the corresponding operator would have dimension 3 x 12; de- 


note this matrix as L 


21 ' 


Now a problem of reduced size can be formulated. The equation 

• ^ 


is constructed, where the correspond to the reduced number of tensor 
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product terms stacked In x^. It will not always be the case that 
Ljj^; in fact, ■ Lj^q and - Lqj^. But the new notation will be em- 

ployed for consistency. The bracketed term In the above equation has dimen- 
sion n X p, with n being the number of Independent state variables and 

a 

p being dependent upon the number of tensor terms retained In the approxi- 
mation. Therefore, x has dimension n x 1 and x^ had dimension p x 1. 
The composition of the vector x^ Is a crucial consideration In the equation; 
a logical scheme Is necessary for consistency, presented In the following. 

Let S be the space of states so that the n-vector x Is an element 
of S. Let U be the space of Inputs so that the m-vector u Is an ele- 
ment of U. Furthermore, let {a, ,a»,...,a } be a basis for the space S 

1 z n 

and (b, ,b^,...,b } be a basis for the space U. The space 
1 z m 

q times 

would have a basis of q-vectors, n*^ in number, of the form 

1 2 q 

where are Integers between 1 and n. Such a space would have 

linear combinations of products, 

12 _ q 

X ® X 9 . . . dl) X , 

of q n-vectors, «:dch comprising a n*^-dlmensional object. But, as seen in 
the discussion above, this number can be reduced If redundant terms within 
the product are eliminated. The number of distinct elements In the product 
Is given by 

,n + q - 1. 

^ q ^ ’ 
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that is, "The combination of (n + q - 1) items taken q at a time." Like- 
wise, to span the space 


U 9 U • ... ® U 

r times 


If * 

m r-vectors of the form 


b . ® b . ® . . . ® b , , 

3l J 2 


where j 2 » • • • » integers between 1 and m, would be necessary. But 
the number of distinct elements in a product 


would be given by 


u^ ® u^ ® ... ® u^ 


,m + r - 1. 

( ^ ). 


The extension of this to product spaces of the form 

^ 0 s ® ... ® s ® u ®^. . . ® u 

q times r times 

follows naturally. The number of (q+r)-vectors of the form 

a , ® . . . ® a . ® b . ® ... ® b , 




1 


‘1 *q -'1 

in the basis is given by the product of n*^ and m . Here, the number of 
distinct elements in a product 

1 <1 1 A 

X ^...®x ®u ®... ® u 

is calculated according to 

,n + q - 1. ^m + r - 1. 

< , ’ • < r > • 

Knowing the number of elements in Xj^ is half the problem; some or- 
dering convention need be established. To settle this issue, consider a 
simple illustration with n ■ 3 and m - 2; thus, {aj^,a 2 ,a 2 > forms a 


basis for S and {b^,b 2 ) forms a basis for U. Then» for the product space 


S 9 S (S ^ f 

a legitimate basis would be 

{a^ ® ® ^1* ^1 ® ^1 ® ^2* ^1 ^ ^2 ® ^1* ^1 ® ^2 ® ^2* 

4^ O ® ® ® ^2* ^2 ® ® ^X* ^2 ® ^X ® ^2* 

*2 ® *2 ® '’x' *2 ® *2 ® *’2’ *2 ® *3 ® '’x* ®2 ® *3 ® '* 2 ’ 

® «i ® 1>2^» *3 ® *1 ® b,, a, ® a, ® , a, ® a, ® b,. 


2’ “3 


'X’ "3 


43 ® a^ ® bj^, a^ ® a^ ® b^} , 


2 X 

consisting of X8, or (3 ) * (2 ), members. The space const' is of Xinear com- 
binations of products of the form 

X ® X ® u 

which wouXd aXso consist of X8 elements; it is the ordering of these ele- 
ments which is the primary concern here. 


Note the manner in which the members of the above-mentioned basis are 

listed. By this convention, the index of each item in the general form 

a ® a ® b 
ij_ 1-2 h 


is Incremented through, beginning at the right. For Instance, with 1^ ■ X, 
I 2 * X Initially, is Incremented from X to r, at which time the in- 
dex I 2 is incremented by one and is Incremented from X to r again. 

This process is continued for i 2 ~ X,2,...,q, at which time 1^ incre- 
ments by one and the previous steps are repeated. Thus, with ij^ ■ X,2,...,q, 

the last basis element in the list will always be a ® a ® b . The 

q q r 

(n^) * (m^) elements of the tensor product x ® x ® u can be listed in a 
similar fashion, the obvious difference being that these elements are or- 
dinary products of the members of x and u. To eliminate redundant terms 
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In the tensor products for identification of the corresponding requires 

a slight modification of this scheme. 


The algorithm for ordering the non-redundant terms from the tensor pro- 
ducts will now be presented via a flow diagram, represented in Figure 2.2. 
Variable names for this flow chart are detailed in the following: 

X and U: arrays of N states and M Inputs, respectively. 

INDX: array of Indexes for elements of x, ordered. 

INDU: array of Indexes for elements of u, ordered. 

NXTR(l): (marker) cumulative number of elements in INDX 

for tensor products of x up to the 1-th degree. 

NUTR(l): (marker) cumulative number of elements in INDU 

for tensor products of u up to the i-th degree. 

XPROD: array of products of combinations of elements of x. 

UPROD: array of products of combinations of elements of u. 

XPR(l) : (marker) cumulative number of elements in XPROD 

for tensor products of x up to the 1-th degree. 

UPR(i) : (marker) cumulative number of elements in UPROD 

for tensor products of u up to the 1-th degree. 

XL: array containing all non-redundant product terms of 

states with Inputs, ordered. 


Three basic divisions comprise the algorithm. The first, represented 
in Figure 2.2a, establishes the order of the Indexes or subscripts of x and 
u to be used in forming products. This first stage is shown in a long form 
for clarity, since a general, condensed version would be less illustrative; 
certainly one set of nested loops would serve the general purpose as op- 
posed to a set for each degree of tensor products as shown. To exemplify 
the ordering of these indexes, consider the previous case of the 3-vector 
X. For an approximation including third degree tensor product terms the 
array of indexes for x, ordered according to the left branch of Figure 
2.2a, would be given by 

INDX -d,2,3;ll,12,13,22,23,33;lll, 
112,113,122,123,133,222,223, 
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2 3 3 , 3 3 3). 


In the above, the notation is such that a semi-colon separates groups of in- 
dexes corresponding to the different degree of tensor product, and a comma 
separates those Indexes which are to be used in ordering product combinations 

a 

of x^f X 2 # *nd Xy Thus, stage two (Figure 2e2b) uses these indexes to 
compute combinatorial product terms; the result for this case would be 

AifROD * ^^1 * ^2 ^ ^3 * ^ 1^1 * ^ 1^2 * ^ 1^3 ^ ^ 2^2 * ^ 2 ^^ * ^ 3^3 * 
x^x^x^ , x^x^x^ , * ^ 1 ^ 2^2 * ^ 1 ^ 2^3 * 

x^x^x^ , x^x^x^ , X2X2X2 , x^x^x^ , x^x^x^) . 


Similar calculations are performed for the m-vector u (as presented, the 
algorithm assumes m > 0), so that in the final stage. Figure 2.2c, the 
product terms are combined and ordered in the vector x^, corresponding to 
the order chosen. 

These final two stages of the algorithm, as represented in the flow- 
chart, are given for the general case (arbitrary degree of approximation). 
The construction is such that the number of multiplications performed is 
held to a minimum. Furthermore, the division of the algorithm into three 
stages offers some conceptual convenience as well. It is not difficult to 
see from this a method for implementation on a digital computer. 


The heart of the identification scheme, then, is embodied in the 
cons tij-t ion of the p-vector The following example will serve to 
summarize the principle. 
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Consider the nonhonogeneous vector differentiel equation 

X • f(x,u), 


with some initial condition x(0) • 
given by , - 



Suppose that in the approximation up 
terms are kept. Thus, the expansion 
given by 


:q. Let the state and input vectors be 



to (and including) third degree tensor 
with tensor notation applied, is 


f(x,u) - 4o ^01 “ 

+ L 2 Q X ® X + X ® u 

+ I<02 “ • “ ^30 * • * 

L 22 X • X • u L^2 X • u • <1 

+ Lq2 u • u • u . 


Putting this into the form discussed above gives 

^01 4o 4l ^02 ^30 ^21 ^12 ^ 03 ^ \ * 

The number of elements in is calculated according to 



where each p^ 


corresponds to one of the linear operators 

p - * ('r') * 


Then, 
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, ( 3 * 1 - 1 ) . ( 2 * 1 - 1 ) , ( 2 * 2 - 1 ) 


• Z + 2 ■¥ 6 


+ 10 + 6 


+ 3 * 3 + 4 


Rathtr ch«n list all of cha 55 mtmbars hara, conaidar Juat a ftw of ':ha 


partltiona co ba atackad in tha vac tor x^: 


iron 



from 


*l’'l“l 

*l*l“2 

*l*2“l 

*1*2“2 

*l*3“l 

*1*3“2 

*2*2“l 

*2*2“2 

X2X3U1 

x.x.u. 
2 3 2 

x.x.u. 


from 


*l“l“l 

*lV2 

*l'*2'*2 

*2“l“l 

*2V2 

*2'*2'‘2 

*3“l“l 

X3U1U2 

*3“2“2 



This illustration should furthsr sstabllsh ths ordsrlng convantion. 


With tha construction of ths aquation complats, ths buslnass of idan- 
tiflcation of tha nay ba undartakan. Tha n sinultanaous nonlinaar 

diffarantial aquations ara intagratad and tha rasulting data is saaplad at 
tha salactad h tins points. Thasa samr.lc^ valuas ara loadad into tha 
matrix which is now pxh in sics. Nota that tha first n*HB rows 

of ara datamlnad from tha samplad valuas of x and u; tha ra> 
malning p>(n4m) rows ars nultlplas and combinations of thosa first n4m 
rows. Fiiully, tha X matrix • dimansion nxh, is formad by loading dar* 
ivatlva astlmatas for x. • x..,...,x at tha h tins points. 

1 A A 

With this, ths matrix aquation assumas tha form 


* ■ ^^10 ^01 ^20 ^11 ^02 ^30 ^21 ^12 ^ 03 ^ \ ’ 


I 


pxh 


nxh nxp 

with matrix dinanslons at shown. Tha laast squaras minimisation algorithm 
is again usad In Idantlfylng tha nxp partitlonad matrix containing tha 
daslrad Tha matrix raturnad will hava flva partitions If only 

sacond dagraa tarms sra kapt, nina partitions If third dagraa tarms ara 
addad, and so on. Tha sorting of the partitions Is martly a numerical 
bookkeeping Job, since thalr sizes arc known in advance. 


As noted at tha beginning of this section, more time points ara re- 
quired for sampling In the Identification now, due to tha behavior of f.hc 
forcing functions. Whereas previously 15 sampling times ware usad, not 
all evenly spaced, now h«40 points arc usad with a constant sampling 
rata. Tha number of points and the rata may ba varied according to the 
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d«BMods of Cho probloa. 


Tho next ooction lllutcracos ch« application of cha aofevaro on aavaral 
raproaoncatlva oxaaplaa, aach ona prograsalng In eoaplaxlty. In chooalng 

cha nonhoBOganaoua vactor dlffarantial aquation for an axaapla, ewe con- 

* 

diciona nuat ba aae: 

1) ause ba aeabla at cha origin, chac la, hava 
aiganvaluaa with oagacl>^a raal parca; 

2) Tha origin muse ba an aqulllbriua point; 

Chat is, 

fj(0,0) - 0 

for 1 ■ 1,2,. . . ,n. 

Cholca of an initial condition is not Intuitivaly obvious in thasa nonllnaar 
axaaplas, so in all casas small partubacions from tha origin ara chosan 
as Initial conditions. 


2.1.3 Second Order Examples 

This section treats two examples, aach of which involves a state vac- 
tor of two alaments and an input vector of two elements. Keeping up to 
third degree tensor terms, the nonhomogan::ous , nonlinear vactor differential 
aquation 

X ■ f(x.u) 


is approximated by 


* • 4o * ^01 “ 

+ L 20 * • X + Lj^l X • u + Lq 2 u « u 
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■f X O X O X + L2^ X o X O u 

+ L^2 X o u 9 u + u o u O u . 

In each case» for a given operating point (x^.u^), the nine individual 
linear operators are identified, and then used to reformulate the two equa- 

a 

tion system to be Integrated for comparison to the true solution of the or- 
iginal equation. 


For the first example, consider the system of two exact nonhomo geneous 
differential equations 

fj^(x,u) « Xj^ 

2.2 

- Vl + "2 - *1 ’ 


■ *2 


where the state vector is 


and the input vector 




X - ( Xj ^, X 2 )' 


u » (uj^,U2) ' . 


The initial conditions to be used in the identification of the L^, are 

jk 

given by 


Xj^(O) » 0-2 , 

X2(0) » 0.4 . 


Notice that the condition of stability is met, that is 
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I has eigenvalues with negative real parts. Furthermore, 

f,(0, 0) - 0 

j 

^ for 1 * 1,2, that is, the origin is an equilibrium point. 

The input forcing functions chosen are sinusoidal in nature, and are 
given by 

Uj^(t) ■ sin(27Tfj^t) , 

U2(t) ■ (0.5) sin(2Trf2t) , 
for f ■ 5 hertz and f 2 ■ 10 hertz. 

A fourth-order Runge-Kutca routine is employed to Integrate the vec- 
tor differential equation; the same routine (and integration stepsize) 
is employed later to integrate the system as embodied in the identified 

model. Thus, a block of data the "true solution" is established and 

stored in arrays corresponding to the and X 2 solutions. The 

stepsize in the integration is taken to be 0.005 so as to comprise 
ample information of function behavior for the derivative estimation. 

Note that for a frequency of 10 hertz for U 2 , corresponding to a period 

.1 
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at *» ' ■»■ *»> ■“ ^ 


of 0.10, each cycle Is visited 20 times In the Integration. 

With the data tabulated, the sampling operation ensues. The true so- 
lution Is sampled at 40 points In time at a constant sampling rate of 40 
samples per second, corresponding to a sample period of 0.025. 

Thus, for the matrix equation 

^ ^01 ^20 ^11 ^02 ^30 ^21 ^12 ^03^ \ ’ 

the matrix has 34 rows and 40 columns. As explained In the previous 
section, the number of tensor product terms (rows of X^) Is calculated 
using the general formula 

n+q-1 m+r-1 

P. - ( ) • ( ) 

q r 

for the q occurrences of x and r occurrences of u In the 1-th 

product. Then the total number of rows of X^ Is found according to 

9 

P ■ I Pi » 

1«1 

and, for this particular example, 

p-2 + 2 + 3+ 4 + 3 + 4 + 6 + 6 + 4 

“ 34. 

The 40 columns of X and X^ correspond, of course, to the number of 
time points used In the sampling. 

The data for the X matrix (dimension 2 x 40) Is obtained directly 
from the (stored) true solution, where the derivatives are estimated ac- 
cording to 

x.(t,+0.005) - x,(t. -0.005) 

fil . -L..J LJl 

^*^lj 2(0.005) 

for 1 ■ 1,2 and J ~ 1,2,..., 40, where the t^ are the points In time 
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chosen for sampling and represents the Ij-th element of X. 

Having completed the sampling, derivative estimation, and loading, the 
least squares minimization algorithm (routine SIMEQUAT) is executed on the 
matrix equation. Knowing the sizes (2xp.) of the individual linear op- 

a ^ 

erators facilitates a partitioning of the 2x34 matrix returned in the 
identification scheme. In Figure 2.3 these operators are shown with the 
corresponding partition of the x^ vector. The eigenvalues of are 

both negative and real. It is interesting to attach some meaning to a 
few of these numbers in relation to the original equations. For example, 
note the relative size of the (2,2) element of corresponding to 

the presence of term x^X 2 in the second equation; the relative magnitude 

<\» 2 
of the (1,3) element of indicates the occurrence of the U 2 term 

of the first equation; Che relative magnitude of the (2,2) element of 

2 

represents the occurrence of the term u^U 2 in the second equation. 

The next step is to integrate the system with the coefficients as 
Identified in the third degree approximation. To do this, the equation 

* “ ^01 4o Hi hz Ho Hi H 2 H3^ *l 

is reconstructed so that the system of two differential equations, each 
with 34 terms in the sum, may be solved yielding x. In the first analy- 
sis, the same initial conditions used in the identification of the iL.. 

jk 

are employed to simulate the original system, with the sane forcing func- 
tions as used in the Identification. The results of this test are de- 
picted in Figures 2.4 and 2.5. As will be the case in all plots to fol- 
low, curve A represents the "true solution" of the original system as 
given by the fourth-order Runge-Kutta algorithm; curve B represents the 
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simulated solution which again employs the Runge-Kutta routine to integrate 
the system as embodied in the identified Figure 2.4 shows the tran- 

sient region of solution x^; the effects of the input functions are evi- 
denced by the oscillations. The simulated solution tracks the true solu- 

s 

tion quite well in this region. Figure 2.5 shows a plot of X 2 from time 
0.0 to time 0.5. Again, curve B overlays curve A. Since the model sys- 
tem was Identified from samples taken out to one second (i.e., 40 time 
points spaced at 0.025), an obvious question would concern the performance 
of the model system for simulations beyond the transient region, out to 
four seconds. The plots for these simulations are given in Figures 2.6 
and 2.7 for and X 2 » respectively. These plots indicate that the 

model system simulates the original system well, far beyond the interval 
in which samples were taken in the identification procedure. As might be 
expected by Inspection of the original equations, after transients have 
settled the forcing functions tend to dominate the solution; this is de- 
picted by the oscillatory motion about the origin for X 2 in Figure 2.7. 

At this point it is interesting to discuss the concept of an "oper- 
ating region", or feasibility region for initial conditions of the states 
of the system. For Instance, the initial conditions chosen for this ex- 
ample were 

x^CO) - 0.2, 

X2(0) - 0.4, 

and the nine linear operators for the series expansion were identified ac- 
cordingly. But suppose that the model system as identified for the above- 
mentioned luitlal condition were used to simulate the original system at a 
different initial condition, namely 
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Figure 2.4 



Figure 2.5 





OWG|f<4L 
Of POOR 


fAGB IS 

WAirv 


Figure 2.6 



Figure 2.7 


x^(0) - n.3, 


X2(0) .’,3. 

Th«t« Initial conditiona repraaant a ehanga in magnituda of SO and 25 par- 
cant, raapactivaly, ovar tha original initial conditiona. Tha raaulta of 

a 

thia taat ara givan in Figuraa 2.8 and 2.9 for and X 2 . A compariaon 
of Figura 2.6 with Figura 2.8 ahova a alight daclina in tracking accuracy 
in tha lattar for tha incraaaad initial condition. But tha plot for X 2 
axhiblta aaaantially idantlcal accuracy aa tha almulation for tha ayatam 
aa Idantlfiad, ainca tha Initial condition waa dacraaaad in magnituda. 

Bafora proceading to dlacuaalon of tha naxt cxampla, an obaarvatlon 

concaming modal aiaa la in order. Sinca thara ara only two tarma of da- 

2 2 

graa higher than two in the ayatem, x^u^ in the flrat quaation and u^U 2 
in the aacond, one would auapect that an approximation keeping aacond da- 
gree tenaor terma would aufflce. Slmulatlona for that Identified ayatam 
proved to be aa accurate aa the former approximation which kept up to 
third degree terma. Since a second degree approximation Involves con- 
siderably less calculation (five partitions, 14 members in x^), one 
might wish to consider using it in this case. 

The second example in thia section is more complex than the first in 
Che sense chat the nonhomogeneous system of differential equations is not 
exact, but a combination of hyperbolic functions. Again, intuition fails 
somewhat as to the behavior of the nonlinear functions, particularly in 
regard to choice of Initial conditions. 

Consider the following system: 

Xj^ - fj^(x,u) 
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Figure 2.8 
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Figure 2.9 


2 "l 

■ u^U2 coth(X2X2) * • «lnh(2x2) ~ 5 sinhx2; 

*2 ■ f2(x.u) 

“l“2 2 

■ • finhCx^) - U2 co«h(x2X2) + slnh(x2)« 

Initial conditions used in the idtntificstion will bs 

Xj^CO) - 0.05, 

X2(0) - 0.075. 


The forcing functions ussd in ths initial integration avs given by 

u^Ct) ■ 0.2 sin(2irf2t), 

U2<t) ■ 0.3 Sin(2irf2t). 


where f^ 

• 5 hertz and 

h 

- 8 

hertz. 

Note that the tt>o conditions de- 

mended of 

the example are 

satisfied; 





f(0. 

0) - 0 

$ 

and 










3X2 



ho 


®^2 

Hi 













x • 0 






u ■ 0 





-A 

! 


‘J 

has eigenvalues in the left half complex plane, thus satisfying the sta* 
bility requirement. 

To begin the identification scheme, the original system is first ln~ 
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tegrated, with a stepsize of 0.005 again, given the above-mentioned operating 
point and Input specifications. From this block of data, 40 samples are 
taken at 0.025 Intervals, and third degree tensor terms are retained. In 
this case, the partition returned by the least-squares identification, 

y ,-1.526 -8.210, 

10 ^ 1.005 1.262^ 

is similar to the amlvtlcal result derived above, and has eigenvalues with 
negative real parts. Vhen the identified system is sir i! ited (again, it con- 
sists of two eqxiations with 34 terms each in the sum) , the results are very 
good for the region up to one second, as depicted in Figure 2.10 for 
and Figure 2.11 for This behavior is expected since the true solution 

was sampled on exactly this Interval. But an extended look at these solu- 
tions shows an unacceptable tracking error beyond the transient region. 

In fact, for this system (as compared to the first example) the transient 
behavior takes longer to settle out due to the relative size of the elgen- 
values of 

The next step in the overall identification, then, is to Increase the 
interval over which samples are taken so as to encompass complete informa- 
tion of the transient region. This can be done in one of two ways; either 
the sampling rate may be decreased (sampling period Increased); or, the 
number of samples taken may be increased. Accuracy might suffer with the 
former, while the latter choice could Increase computational complexity. 

When the number of sample points is Increased to 80 with the same 
sampling period of 0.025, the region of the true solution to be observed 
is from t *> 0 to t ■ 2 seconds. Thu.'i a larger portion, if not all, 
of the transient region would be sampled. Using the same initial condl- 
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tions and the same specifications for the control inputs, the first parti- 
tion of the identified system 

2. f-2.023 -5.093, 

Ho " ^ 1.009 0.941^ 

shows much ‘better agreement with the analytical expression for L^^q. Fur- 
thermore, the tracking is greatly Improved, as Illustrated in Figures 2.12- 
2.13 for and ^ 2 * respectively. Note that in each solution curve B 
overlays curve A exactly for the region up to two seconds. But there is 
a slight tracking error beyond this point which remains approximately con- 
stant in time after three seconds. This constant steady-state error would 
be acceptable in many applications, although more accuracy is possible. 

A sensitivity analysis on this system as identified for Xj^(O) ■ 0.05 
and X2(0) ■ 0.075 shows good results when the initial conditions are 
varied In magnitude, that Is, 

Xj^(O) - 0.10, 

X2(0) ■ 0.03. 

The results using this initial condition (with the model system identified 
for the original initial condition) to simulate a corresponding true solu- 
tion are given in Figures 2.14-2.17. For the variable Xj^, Figure 2.14 
shows a good overlay of curve B onto curve A, with a small error 
arising around the peak in the solution. Figure 2.16 shows a plot over 
a wider range of time, and reveals that the simulated solution (for x^ 
again) actually Improves past four seconds. Solution X 2 * represented 
in Figures 2.15 and 2.17, does not exhibit quite as good results, but still 
maintains an acceptable curve fit with constant error. 

The accuracy of these simulations can be Improved upon even more by 
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sampling the original system over a yet wider region. So consider a sam- 
pling scheme which again uses 80 sample points, but which now takes sam- 
ples at a rate of 25 per second as opposed to 40 per second. This cor- 
responds to a sampling period of 0.040 seconds so that for 80 samples the 

s 

region from t*0.0 to t~3.2 seconds is sampled. Figures 2.18 and 
2.19 show solutions and X 2 out to four seconds for the model system 

using those initial conditions as used in the Identification (again taken 
to be Xj^(O) ■ 0.05 and X2(0) • 0.075). The tracking is very good for 
this system now, as in each solution curve B overlays curve A through- 
out. Now a sensitivity analysis is performed by simulating with the system 
model using x^CO) • 0.10 and X2(0) « 0.03. A comparison of Figure 2.16 
to Figure 2.20 for solution reveals that for the latter plot, which 
corresponds to the model system identified from samples taken to 3.2 sec- 
onds, the match of the simulated solution to the true solution is much better. 
This is partictilarly noticeable in Figure 2.21 for X 2 , upon comparison 
to Figure 2.17 for the previous identification. It would seem that even 
more accuracy could be achieved, although this last pair of plots indicates 
that a good system model has been achieved without suffering any additional 
ill effects in computational complexity. In fact, it is evident that the 
decrease in sampling rate in the last identification caused no noticeable 
error; to the crntrary, it led to more accurate results. 

These two second order examples presented in this section were dis- 
cussed in detail to Illustrate the flexibility of the method. The first, 
an exact vector differential equation, was chosen to give some meaning to 
the elements of the individual t.. operators and to Illustrate that the 
degree of the approximation is in general problem intensive. The second 
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Figure 2.12 
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Figure 2.21 


«xampl« shoved that even a very unattractive nonlinear, nonhomogeneoua 
system can be adequately handled by this technique. 


2.1.4 Third Order Example 

This qectlon treats an example given by a system of three nonhomo- 
geneous differential equations. The state vector. 


X ■ (Xj^, Xj, Xj)*.. 


consists of three elements, as does 


u(t) - (uj^(t), U 2 (t), UjCt))*, 


the input vector of forcing functions. Obviously, the size of the problem 
will be Increased compared to that of the second order systems of two states 
and two Inputs. Whereas for n ■ m » 2 In the previous examples the length 
of the tensor term vector x^ was 34, now the length of x^^ will be p • 83 
for an approximation keeping up to third degree terms. Again, this number 
Is calculated according to 


P - I 


9 

1-1 ^ 

for the 1 partitions. Expanding this for n * 3 (three states) and 
m - 3 (three Inputs), the result Is 

p-3 + 3 + 6 + 9 + 6 + 10 + 18 + 18 + 10 
- 83. 

Thus, the matrix equation to be constructed takes the following form with 
matrix dimensions as Indicated: 


^ “ ^^10 ^01 ^20 ^11 ^02 ^30 ^21 ^12 ^ 03 ^ \ 

I I V 

3xh 3x83 83xh 


Here, h Is the number of sample points used In the loading scheme for 
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estimation of derivatives and construction of X^. 

The third order system of this example Is of the same nature as the 
previous example In that exponential and hyperbolic terms comprise the non- 
linear system. Consider the following equations, 

» 

fj^(x.u) - Xj^ 

u U- 

• Uj^e coshCxj^Xj) - e slnh(3Xj^); 

f^Cx.u) - X2 

2 2 

■ Uj^ cosh(xj^) - 2 slnh(x 2 ) + Uj^u^ cosh(x^) ; 
f^Cx.u) ■ Xj 

■ Uj^ slnh(xj^X2) - slnhCx^) . 

To Identify die model system, the Initial conditions chosen are given by 

x^(0) - 0.05, 

X2(0) - -0.05, 

X3(0) - 0.08. 


Choices for forcing function Inputs are sinusoids with the specifications 

Uj^(t) ■ sln(2irfj^t), 

U2(t) ■ (0.5) sin(2Trf2t), 
u^Ct) • (0.5) sln(2irf3t). 


where the frequencies for each Input are given by 


In hertz. 


- 1 , 

- 5 , 

- 10, 
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This system Is constructed such that the following conditions ere 


easily verified: 


f(0, 0) - 0 ; 


and for stability 




" r — ^1 
i 


X 

u 


0 

0 


-3 


0 


0 

0 


-2 


0 


0 

0 

1 


With the problem thus formulated the technique Is applied. Using an 
integration stepsize of 0.005, the system Is integrated and the true solu- 
tion stored. This data Is sampled at 40 time points evenly spaced at In- 
tervals of 0.025 for loading of the tensor term matrix and for derivative 
estimation. Thus, the sampling takes place over the first one second In- 
terval of the solution. The first partition returned by routine SIMEQUAT 
In the least-squares minimization is given by 



-3.072 -0.112 0.440* 

0.009 -1.948 0.003 

-0.001 0.005 -0.936 


which has eigenvalues with negative real parts. Note the similarity to 
the analytical expression for given above. 

Simulations using this model system consist of an integration of the 
three equations which have 83 terms each in the sum. Over the Interval in 
which the 40 samples were taken, the simulated solution matches the true 
solution very well when the initial conditions are those as used in the 
Identification. These results are shown in Figures 2.22-2.24 for x^,X 2 > 
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and raapactlvaly. Flguras 2.25>2.27 depict tha aaoc ainulationa 

over a wider range of time. The most error occurs at the peaks of the so- 
lution curves, especially noticeable in the variable (Figure 2.25). 

Consider now a sensitivity analysis on this model system for the 
following changes in the initial state conditions: 

Xj^(O) - 0.10 , 

XjCO) - -0.10 , 

XjCO) - 0.16 . 

Using the same forcing function as in the identification, the original sys- 
tem is integrated with these new initial conditions to form the true solu- 
tion. Then the model system is employed with these Initial conditions to 
simulate the true solution. These results are shown in Figures 2.28-2.30 
over an eight second Interval of time, well beyond the transient region of 
each solution. Very little error results between the simulated solution 
and true solution. 

A conclusion of the final three plots is that for this particular ex- 
ample an ample number of sample points for system identification is 40; 
perhaps even fewer would produce equally acceptable results. In the final 
analysis the Initial condition of each variable is double that which was 
used to Identify the model system, indicating (by the accuracy of the 
tracking in the simulations) that an even greater variance in the initial 
conditions would produce acceptable results. Thus, this model system seems 
to have a larger feasible operating region than that of Che examples treated 
in the previous section. 
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Figure 2.22 
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2.1.5 Conclusion 


An important aspect of the identification scheme as put forth in this 
chapter concerns the capabilities of routine SIMEQUAT and the least squares 
technique in general. A question of conditioning arises as the size of the 
problem increa&es. Three factors contribute to problem size: the length of 

the state and input vectors; the degree of the approximation; and the number 
of sample points used in observing the system* For larger problems it is 
this last factor which in some sense limits the actual identification algorithm. 

Some comments concerning actual computer time are in order. As an ex- 
ample, the identification for the third order example presented earlier was 
carried out using interactive software on a terminal over a time-sharing net- 
work, IBM system 370. The entire identification process took approximately 
ten minutes terminal time, corresponding to about five seconds of CPU time. 

The verifying simulations took almost as long, due to the size of the re- 
constructed equations. Smaller examples are even faster, although the 
above-mentioned times are certainly reasonable. 

As a firial note, it is interesting to consider use of other sampling 
techniques in the identification procedure. For instance, suppose that 
variable sample periods were used in an attempt to better observe the sys- 
tem. In fact, an apriori knowledge of particular behavior of the observed 
data could lead io clever cho ^es for sample periods. The time between 
samples could even be generated randomly. 

Results of the discussion and examples of this chapter indicate that 
the technique can be utilized in a number of applications. Given a block 
of observed data from a nonlinear system, and given knowledge of the num- 
ber of states and inputs, several model systems (for different operating 
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regions) in the fora of e collection of Indlvldusl llnesr operators can be 
Identified. 
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2.Z Report on CARDIAD Progress (R.M. Schafer) 

'i'toug’i funded work on the CAROIAD method was discontinued more than 
six months Its progress has continued steadily by means of fellowship 
support extended by the University of Notre Dame. During this period, the 

t 

emphasis has been placed upon developing the method to be helpful with di- 
agonal dominance designs in cases having four or more inputs and outputs. 

Appendix B contains the preprint of an especially challenging design 
on such a model. 

2.3 Report on Feedback Loop Closures (V. Seshadrl) 

Also mentioned in the last Semi-Annual Status Report was a recent ex- 
tension of early grant work in the general area of pole and zero assign- 
ment and the exterior algebra. This work received support from Grant 3048 
in its early stages, but in later years received assistance from different 
sources . 

During the period of this report, another presentation has been made 
on this subject. The preprint is contained in Appendix C. 
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cuPiAP AFPiOAOi TO STSTZM Mcnuna 
vm AriLxcAixoi to tdibofan Bfcmt msdkls* 

X. mehMl SctaAftr aad KichR«l K. S«la 
D«p«rtMat of Uoetrieol Titgln— ring 
DblToroity of llotro Ooat 
Koero Dmo» TndlaiM 
m 46536 



Om of tbo f ootoroo of prooonc dmy rooooreh on 
liaoor ■ulcivoriabXo wymtmM boo boon o roumod laeor- 
•ot la Kpqttlot Mtbodo. laood ttpoa tbo dotoalaaac of 
rotitra dlff ttoaeo, thooo aocbodo aaoc dovolop proeo* 
dtiroo which iatorfoco with star oyMOCric, aultllixMar 
foru. ▲ woU teown latorfoco hoo boon aodo bp loooo- 
brock, who uood a coocopt callod diagonal doalaanea. 
This papor raporta on a graphical, latarmetiTo wap to 
achiowa cha concapc. 


Introdnetion 

Za meant paara, laeraaaad attaatlon haa haaa 
paid to tha uoa of fraquanep doaaia tacfaoiqaaa for tha 
daaign of ■altiwaciabla control apataaa* Moot of 
thaaa taefaniquaa ara b a a ad upon tha aqnaclon 

p^(a) - |M(a>| PqCo) (1) 

which ralacaa th^ aaroa of tha opon loop charactaria* 
tie polpnowial p^Ca) aad tha laroa of tha cloaad 
loop charaetariat^ polpnoalal p^Ca) through tha da* 
taralnant of tha ratuxn diffaranca aatrlx M(a). Ciran 
that tha aaroa of tha opan loop charaetariatlc polpnow- 
lal ara taiown, acabilitp of tha cloaad loop charactar* 
iatic polpnonial can ba dataralaad bp Vpquiat analpoia 
of |M(a)U U&fortunatalp, dlraet Hpquiat analpaia of 
|M(a) I plalda littla daaign iMight. Tharafora, al* 
taraata naana of acndplAg |M(a) | hawa boon dariaad. 

In tha Isvaraa Npquiat Arrap approach dua to 
Koaanbreck, tha apatm la firat eoapanaatad to aehiawa 
diagonal doalaanca. An n z n aatrlx Z(a) la aald to 
ba dlagooallp coltan dpi inant if for all s on tha 
Npquiat contour 0, aad 1 • 1,2, •••,n, 

!««(•) I > I W 

“ j-i 
iH 

If thia eondltlon la aaclaf lad, tha uaual not anelrcla* 
nanta nada bp ba Npquiat plot of |M(a) | ara aqual to 
tha iun of cha nac anclrelananta nada bp Cha diagonal 
ontriaa of M(a). Thu, atabilltp can ba datamlnad bp 
Npquiat analpaia of tha diagonal ontriaa of M(a). 

Tha CAtDIAS (Coaplas ^eoptabiliCp kagion for OIA* 
gonal naninnea) aatbod la*a graphical cachnlqua for 
aehlv^ng chla daninnea condition. 

CAPIAD Wathod 

Conaidar tha ajmtm of Plgura 1. Tor tha purpoaaa 
of thin papar, 6(a) rapraaanta a 4-lapuc, 4««otpot 
nodal of a turbofan jot angina. It la daalrad to da- 
aign tha eonpanaator K(a) aoeh that 6(a) K(a) la 
eoluan doniaant. Tha eonpanaator la aotnallaad to 
hawing I* a on cha aaln dl^onal, ao that donlnanea la 
aehlvrad la a glwan colnan of 6(a) C(a) bp approprlata 


*Thla work has baan aupportad in part bp cha National 
Aarooaueica and Spaca idnlnl aeration undar 6rant HS6* 
3045. 


eboiea of tha off-diagonal antrlaa in tha corraapon* 
ding ealiM of K(a). 

At a fraquanep, a aufflciant condition for doni* 
nanca can ba onpraaaad in a quadratic inaqualitp of tha 
fom 

f (*> - x^Ar ♦ xS -f c > 0, (3) 

whara A, b, aad e ara raapaetiwalp a natriz, a 
waetor, aad a ofalar foznad bp awaluation of tha plant 
tranafar fuaetioa natrlx at cha fraquonep baiag atudlad, 
whara x la a waecor of tha raal aad inaglnarp parts 
of tha off-diagonal ontriaa of a eol»B of tha e o np a n 
aator, aad wham anparseript t danotaa traaapoaa. 
Soniaanca la achlaw^ bp choosing x aoeh that f(x) 
la poaitiwa. 

Sawaral approachaa ara uaad to chooaa x such 
that f(x) > 0. Slnca it is daalrabla to achiawa donl- 
nanca with as ainpla a eonpanaator as possihla, cha 
gradient of f(z) is takan with raspaet to aach antrp 
x^ aaanning all othar antrlaa ara xaro. lara, x^ 
nap ba undaratood as a pair (r^,l4) conaiatl^ of tha 
raal and inaglnarp parts of aona off-diagonal eonp»- 
aator ontrp. Thia approach, mfarrad to as tppa 1 an* 
alpaia, attaapta to achiawa doniaanea la a eolonn hp 
uaiag oalp oaa aoaaaro, off-diagonal aacrp for tha 
colunn of tha eonpanaator. In tha awont that It la 
inpoaaihla to achiawa donlnanea with oalp ona nonaare. 
off-dlagoaal antrp, tha gmdiaac of f(z) with ra- 
space to all wariabloa la takan. Thia approach la ra- 
farrad to as tppa 2 analpaia aad utilisaa all off-diag- 
onal antrlaa of tha eonpanaator to achiawa donlnanea In 
a eolian. A third naana of choosing cha waecor x io 
uaad in tha mwmt that cha haaalan in tha tppa 2 anal- 
paia approach la indafinica. It la Inown chat ana ao- 
Intlon to naking 6(a) doninant la to eonpan as to with 
cha Inwaraa apstm. Thna, a solution for tha waetor 
X la to chooaa tha waluaa of tha inwaraa apstan at 
that fraquanep, nomallaad to 1 on eba dlagou3. so aa 
to fit cha fom of tha eonpanaator K(a). In tha caao 
whara tha haaalan la nagaclwa daflnlta, this inwaraa 
apatan analpaia, Vnown aa tppa 4 analpaia, pradleta 
cha aana solution aa tha tp^ 2 analpaia plots. 

Tha CAPIAD plot la a graphical rapraaancaclon of 
eba raaulta of tha gradlant analpaia. Conaidar tppa 1 
analpaia of a glwan colnnn. f (0, .•«0,x^,0..) la a 
paraboloid in 3-apaea, and tha walua found hp tha gra- 
diaac analpaia can ba a poaitiwa naiimm, a nagaclwa 
naxinnn, a poaitiwa nlninan, or a nagaclwa nlninon. Za 
tha pesidwa narlw caaa, anp walua of which llaa 
l a al d a cha Intaraaction of f(z) and tha eonplax plana 
z^ will naka f(..0,z^,0, «.) poaitiwa; and doainanca 
will ba aehiarad at tha fraquanep baiag atudlad. In 
tha CAIDTAD plot, thia is rapraaantad bp a solid clrela 
which la tha aolntion of f(..0,z^,0..) ■ 0, aad a '4>* 
at tha wnlna of x^ whara tha gradlant waniabaa, which 
is at tha emtar of tha clrela. In cha caaa of a aag- 
atlwa nlnlann, all walnaa of X 4 Iplng outaida tha 
clrela f(.,«0,zi,0..) • 0 will naka f (•*0,z^,0...) 
poaitiwa. In this caaa, aa 'x' la drawn at cha walua 
whara tha gradlant waniabaa and a daabad eizcla at 
f(...0,z^,0.«) *0. In tha nagaclwa naxi— caaa, ao 
walua of will aehiowa doniaanea; aad a 'A* la 
drawn. In tha poaitiwa nintniin caaa, anp walua of z^ 
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vlU «chl«ir« tiwlBMmm; la tha aoliaai at thla fra^aa* 
cjf aad a *Q* la inm* 

Xa tjy a a 2 aa4 4 aaal 7 «aa« tba aaatar m 

4 r a a a at tha gradlaat valaaa« bat tha eaatar typa aad 
elroia typa ara daeldad hf aahiac a aarat eaaa davla* 
tloa item tha gradlaat valaaa af all hat oaa af tha aa* 
triaa of a; aad that tha raaalaiag Mry la aaalTiad 
la a faahloa aaalogoaa to tTpa 1 aaalTola. 

A C4I9X4S plat raaalta ahaa thla graphical gra* 
diaat lafaraatiaa la plottad 9wm a raaga of fraqaaa- 
alao. rigotaa 2 aad 3 ara tppical CdtSXdD plat# aad 
alU ba ttaad to daaeslha eaapaaaatar daalga. 

ngara 2 la a typa 1 aaalpala plot ahUh coatalaa 
ealy aolld elxtlaa. Xa thla eaaa, thara «lat eoaataat 
raal aalaaa (r^,0) far ahiah lla laaida all of 
tha aolld elrelaa* laaaa, to achlava doodaaaea la thla 
eoluaa at all fra^oaaeiaa, aay aoch ehoiea of ai will 
aofflea, alaao f<.«0,x^,0*.) aill thm bo pooitlra at 
all fraqnaaeiaa. Xa Flgura 3, thara aalata ao aaeh 
eoaataat raal oalua, hat a aiapla first ordar aatry 
which aa a foaetloa of froquoacy traeaa cha eeotaro of 
tho elxcloo can ha oaad* Tfaoa, if tha pXot ia«> 

dicataa that ao eoaataat raal rmliia vlU aehiava dood* 
aaaca, tha ehapa of tha plot guidaa tha daaignar ia da- 
tasaiag a fraqoaacy dapaadaat aatry. 

fnlr ****rr^" 

The aodal oaad ia tha foUoviag daalga aaaapla la 
takaa froa (1]« Xt la a tlzth ordar, i-^laput, 4-oat* 
pat daoeriptiaa of a tarbofaa aagiaa* 

Aa a first atop la tha daalga procadiira, tha aodal 
oaa eoapaaaatad with tha lawaraa ayotea avalxiatad at 
a ■ 0. Flgoraa 4-g ara tha typo 1 aaalyoia plots of 
tho 4,2 aatry, tha 3,4 aatry, aad tha aatlra first col* 
40. Typo 1 aaalyaia of tha first eelaiB iadleatas 
-bat doaiaaaca eaaaet ha achlaoad uaiag only oaa aoo- 
taro, eff-diagoaal aatry* Tha aaaa oaa trua for tha 
third eoloBA* Flgoraa P-U ara tha typo 4 aaalyaia 
plots for tha first eeloBa; aad Flgoraa 12-14 sra tha 
typo 2 aaalyaia plots for tho third eelaaa. 


Xa hath tha a a c o a d eoltaa aad tha fourth coltaa, 
doalBaaca oaa achlovahla uaiag typo 1 aaalyaia aad eoa- 
ataat eeapaaaatloa aa daacrlhad la cha diacuaaioa of 
Flgura 2. Doaiaaaca oaa aehiaoad la eoliH 2 by 
chooalag cha 4,2 aatry to ba -SBO.S* Kota that this 
▼alas llaa vithia all solid elrelaa aad oucaida 
daahad elrelaa. Xa Ilka aaaaar, cha fourth coIobb oaa 
aada doaiaaat by chooaiag tha 3,4 aatry to bo -.5f. 


Xa tha first aad third eoliHLa, it oaa oaeaaaary 
to fit all thraa off^iagooal sicrlaa of tha c oap a n - 
aator to tha ahapa of tho eoatars of tha typo 4 aad 2 
plots, raspaetioaly. Ia aach eaaa, aacoad ordar eea- 
paaaatioii oaa naeaaaary to fit adaquataly tha ahapaa. 
Tha thraa aacriaa ehoaaa for tha first eoluBB oara 

. (.) , - ,01?^ ^ 

.227l-2a^ .23$a 1 




-.044a* - 
.2271-2a^ 


2.30a 

4> .23ta ^ 1 


> (.) . — . 

.2271-2a^ > .23$a ♦ 5 

Tha third coli^ oaa aada doalaaac olth cha foUowlag 
thraa off-diagoaal ncriaa 

V (.) , 

.lfW-2.* ♦ .M37. ♦ 1 


w . wm? • .t za r i i . . 

.im-2.' ♦ .M37. ♦ 1 

t, ,(.) . - . MHt . 

.im-2s^ ♦ .0S37a ♦ 1 

Xt should ba oacad that, la aach eaaa, tha danoa* 
laator polyaoBlsl of tho eoluBB ia tha saBO, tharahy 
baaping tha ordar of tha raaultiag nnap ansa ter oboU. 

Kith eaBpaaaatlsB aa daacrihad abowa, typo 1 ob- 
alyais oaa rapaatad to oarlfy that doaiaaaca has baaa 
achiawad. Fl^aa 15-lS ara a typo 1 aaalyais plot 
froB aach coloBa. Kotc thet ia crary eaaa, tha plat 
pradicta that aa aceopcabla ai>latisa ia tha arigia. 
Siaca typo 1 aaalyais is draoa aasMlng all other off- 
disgoaal mtximm ara aaro, this lapliaa that tha cal- 
iBMO ara aoo doaiaaat, alaca tha plots predict that 
Idaaticy r nap anaa tloa will aehiava doBiaaaca. 


Tha tools aaaocistad with Kyqoiat aaalyaia of (1) 
aro oftaa helpful la tha aaalyaia aad daalga of ailtl- 
warlshXa ayacoBa la cha fraquaacr doaaia. Oaa of the 
oaya to approach the daciga of |M(a) | ia (1) ia by 
aaaaa of tha doaiaaaca Idaaa of KoooBbrock [2]. This 
paper daacrlhaa a grapuicsl, lataractiva procadura for 
attalaing doaiaaaca. For other aaaaplao, aoa [3]. 
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ABSTRACT 

In this paper we establish the connection between linear multivariable 
feedback loop closures and the induced exterior map, the latter being the 
by-product of exterior algebras over the input and output vector spaces. 
This suggests the concept of a sequence of multivariable zeros which should 
enhance the designer's ability to shape multi-output transients. 

Introduction 


Much interest has been evinced in recent years about the exterior al- 
gebra, a structure especially suited for addressing questions related to 
matrix determinants and inverses. It has been shown that this structure 
can be used to solve problems in pole assignment [1,2], in individual zero 
placement [3] and indeed in a host of areas related to systems and informa- 
tion theory [4]. The present paper may roughly be divided into three parts. 
The first part introduces the exterior algebra. The presentation is ex- 
tremely brief due to limitations of space; for more details the reader is 
referred to Greub [5]. The second part of the paper considers multivari- 
able feedback loop closures, the loops being closed one by one, and studies 
their relation to certain entities called numerators of the 'xth kind, the 
latter arising from a frequency design method used in industry [6]. The 
last part of the paper establishes the connection between numerators of the 
Vth kind and the kth induced exterior map, the latter being a by-product 
of exterior algebras over the input and output vector spaces [7], In this 
way we will see that when considering multivariable feedback loop closures 
it would be helpful to view them in the context of the appropriate induced 
exterior morphism. 

The Exterior Algebra 


Consider an F-vector space V. We can construct an exterior algebra 
[5] AV over V. The bilinear operator introduced by this construction 
is commonly called the exterior product or the "wedge" product * , and 
operates as 


(cti a^ + a^) . a3 


Oi a^ 


2 2 


h ' ®3 ^ “4 * 4 ^ " ®3 ®1 ' ^3 ^ “4 ^1 ' * 4 ’ 


where a^, a^, a^ belong Co the algebra AV, and a^, are 

field elements from F. Furthermore, the operator - is skew-symmetric. 
Now consider a map f ; V -► W. If we construct the exterior algebras 
AV and AW over Che vector spaces V and W respectively, the map f 
induces a morphism over the algebras [7], which is just a sequence of 


*This work has been supported in part by Che National Aeronautics and Space 
Administration under Grant NSG-3048 and in part by the Office of Naval 
Research under contract N00014-79-C-0475. 
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maps over Che kth exterior spaces, as shown in Figure 1. 



t 


AW 




A,V 

I 


♦ 

A^W 


* 

A, W . . . 

k 


* Figure 1. The Induced Exterior Morphism, f'. 

A reasonable question that could be asked at this juncture is what 
relevance Che induced exterior morphism has to multivariable feedback 

control design. The answer is that there already exists a feedback control 
design method in industry [6] chat makes partial use of the induced excer> 
ior morphism structure. It is noc clear, however, that the originators of 
Che method are aware of the structure and it appears chat the design method 
could be extended by making fuller use of the sequence of exterior maps. 

In order to address these issues, let us first consider a problem suggested 
by Hofmann, ^ al . [6] and study the connection between feedback loop clo- 
sures and certain entitles tailed numerators of various kinds, Che latter 
arising when the following problem is solved. 

Problem : Obtain expressions for the arbitrary closed-loop transfer func- 
tion after feeding back one or more outputs (measurements) to one 

or more inputs (controls), for the plant shown in Figure 2. 


rj^(s) 

U^(s) 


y, (s) 

r, (s) 

u-(s) 

N(s) 

y2<s) 


— " — > 

• Up(s)_ 

A(s) 

• >^m(s) _ 






Figure 2. A Feedback Loop Closure Problem. 

Comments : Figure 2 shows a multi-input multi-output plant; it has p in- 

puts and m outputs. The reference, input and output may be compactly 
called r, u and y, and considered as elements of R(s) -vector spaces R, 

0 and Y respectively, R(s) being the field of ratios of polynomials in 

s with real coefficients, the denominator being non-zero. The plant is 
nth order, linear and stationary. The expression for the plant in Figure 
2, that is N(s)/A(s), is considered to have been derived, with the usual 
assumption of zero initial conditions, from an input-output plant descrip- 
tion in the s-domain written as 

A(s) y(s) » 3(s) u(s) , 

where y belongs to the R(s)-vector space of dimension m, [R(s)]°, and 
u belongs to [R(s)]P. In the following discussion we shall drop the ex- 
plicit dependence of the variables on s, for notational convenience. 

The outputs y may be expressed explicitly in terms of th.^ inputs u 

as 
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where A Is Che nth degree characteristic polynomial for the plant. If 
A and B are left coprime, 

A ■ dec A. 

N is represented in numerical calculations by an mxp matrix whose ele- 
ments belong to R[s], chat is, they are polynomials in s with real co- 
efficients^; N is Che adjusted plant transfer function numerator matrix. 

Solution ; In Che following treatment let the polynomial matrices A and 
B be expressed in terms of m-lengch column vectors as 

A - .^2 *■* *m^ ’ B ■ [b^ b^ ••• bp]. 

Consider the case of a single feedback loop, that is, one output y^^ 
is fed back to the jch comparison point through a feedback gain g.^, 3 

shown in Figure 3. ^ 


Figure 3. The Feedback Problem With One Loop Closure. 

The relationship between the output y and the reference r may be ex- 
pressed as 

(A + BG)y - 3r, 

where y and u are given by 

y * ^^1 ••• ^ * ^’^i ••• V ' 

G is a pxm matrix of zeros except for a single non-zero element ^ 

R(s), and is expressed in terms of p-length vectors in the same fashion as 
A and B above by 


G» [0 ... Og^O 0] 


where 




is the 1th 


Note thet is the jth elanent in g^, and that g^ 

coluan of G. We can then express the equation 


Vc * '*1 * V. 


■ b.r, +..+ b.r, +..+ b r . 
11 d d p p 


Note that in the above equation a^, k • 1, . . ,m and b^^, k <■ l,...,p are 
vectors whereas yj^, k ■ l,...,in and r^, k ■ l,...,p are scalars, the 


field being R(s). 

We wish to isolate the closed-loop transfer function this may 

be achieved by taking the exterior product of both sides of the above equa- 
tion with the (m-1) -exterior term 


*1 *c-l “ *c+l 


i-1 


(Si + 


‘i+1 


m 


which is the exterior product of all the vectors on the LHS of the equation 
except for a^, whose coefficient y^ is the output of interest in our 
current discussion. Because of the multilinearity and skew-symmetry of the 
exterior product, all the terms on the LHS of the equation will become zero 
except for one term which includes If, at the same time, all the 

references except the one of interest, r^, are held zero, the equation be- 
comes 


‘1 

a. 


K a. ^ 

c-1 
a 


-- *1-1 - '*1 * 


c-1 " *c+l 


i-1 


(^i + b^gji) 


‘i+1 
a 


- S * Vc 


i+1 


>> a 


m 


Vd- 


Both sides of this equation contain m-exterior products of m-vectors; the 
products, therefore, are determinants and hence just field elements. Also, 
the products are skew-symmetric so that a^ on the LHS and b^ on the RHS 

of the equation may be moved into position between while 

retaining the validity of the aquation. Thus we can get the arbitrary 
closed-loop transfer function with a single loop closure from the 

output y^ to the input u^ through tne feedback element gjj^, as 


V 

" c 


^1 — Vi 


c+1 


‘i-1 


'*1 * - *1*1 


m 


®1 ®c-l 


*c+l *1-1 - <*1* '■j*Jl> ■ *1+1 


.. a 


m 


or 


(ai *• 


+ a^ 


S-l ~ ^d ' Vl 


c-1 


^d - \+l 


m 

a 


i-i^_^ 


^+1 


.. a 


.Sll> 


(a. 


.A a a, 4. . . < 

m 1 


i-1 


i+1 


*.*ji> 


In order to interpret the above expression for the closed-loop trans- 
fer function observe that, in terms of numerical calculations, each 

of the two terms in the numerator and in the denominator of the above equa- 
tion is simply a determinant; thus a total of four determinantal calcula- 
tions is involved here. Determinants of this sort recur in transfer func- 
tion expressions with one or more loops closed; they are called numerators 
of various kinds [6]. The kind of the numerator is dependent on the nature 
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of th« mix b«tween the colunms of A and the columns of B In the parti- 
cular determlnancal expression. For Instance, the above equation contains 
two numerators of the first kind, that is, tvo determinants resulting from 
the m-exterlor product of (m-1) columns of A with one column of B, 
namely 


and 


“d 1 


• a. 


S-1 - ^^d 


c+1 


a 


. a 


1-1 


’j " *1+1 


Similarly, Che equation has one numerator of the second kind, that is, Che 
n-exterloc product of (m-2) columns of A with two columns of B, namely 




7,7i 
"d“j 


*1 *c-i." ^d 


*c+l 


« a 


1-1 




1+1 


m 


Then, under the condition that dec A > A, we can express the closed-loop 
transfer function compactly as 

^c ^c^i 

y u , ® 1 1 UjU . 

_£. . d d 1 


A + g.., N ' 
u, 


Thus, In the one-loop case, that is, with exactly one feedback loop 
(from y^ to through as shown in Figure 3), the expression for 

the arbitrary closed-loop transfer function Involves numerators of the 
first and second kinds. It has been found [8] by using a similar exterior 
algebraic mechanism for manipulating the closed-loop transfer function ex- 
pression that in the k-loop case, chat is, with feedback loops from k 
different outputs to k different comparison points, numerators up to the 
(k+l)th kind are involved. Or, the highest kind of numerator in the ar- 
bitrary closed- loop transfer function expression is intimately connected 
with the number of loop closures. 


The Numerator of the kth Kind and the kth Exterior Map 


The arbitrary numerator of the kth kind is formed, as an extension 
of Che definition of nximerators of the first and second kind, by taking the 
m-exterlor product of (m-k) columns of A with k columns of B, as 





‘i,+l 






^ a 




.. a 


m 


We will show in the following that numerators of the kth kind may be 
understood in the context of the induced kth exterior map for the input 
R(s) -vector space U to the output R(s) -vector space Y. That, and the 
fact that numerators of different kinds arc intimately related to the num- 
ber of loop closures, would lead to the conclusion that there exists a 
strong association between the induced exterior map and loop closures. 

Consider the equation relating the output y to the input u, ex- 
pressed as 

y » A B u. 

The vectors u and y belong to R(s) -vector spaces U and Y of dimen- 
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slon p and m raspactlvaly. Thus v« have the morphism of vector spaces 

a’^ B : U T. 

We have seen earlier that such a morphism of vector spaces induces a mor- 
phism of algebras 

(a"^) * : AU ■*> AY 

as shown below 


*(A“^)" : 1 a"^B j(A‘^B)2 


• 




. . A,Y 
k 


Figure 4. The Exterior Morphism Induced By A* a. 

Of particular interest is the map 

(A-S)- , - A^T 

because we Intend to show how this map is related to numerators of the kth 
kind. In order to try and compute this map, we express [8] 


(A“^)^ - - (A^)S 


-1 Si 

^k det a; * ^k 
Tc 


Let us assume for the moment that we can find a matrix T such that 

T AT - I det A. 

tc 

The reason for this assumption will become clear in the following; but for 
the moment, let us rewrite the above expression as 

Comparing the above with the definition of adj as 

adj AT 

det A^ * ^ ’ 

and because A^ is invertible, we have 

det A^ det A 

Thus we can replace adj A^/det Aj^ by T/det A in order to rewrite 

« , TBC 

The Nature of the Matriic T 


We now investigate the nature of the matrix T; to understand T is 
to understand the nature of the kth exterior map (A**^B)r, because 
(A*^B)^ is just TB,/det A. Recall that T has been defined in terms of 
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its acclon on th« iaag* of A^, as 

T • I dat A. 

Nov A^ is caleulAtad by forming all poaslble lock minors from the matrix 
A, and is thus 

^..k,l-*it •* ^..k,i^..i^ *’ ^..k,m=^..m 

I 


|j‘m-kfl..m»l..k ** "m-fc+1. .m,i^. .i^ " "m-fcfl. .m,n-k+l. 

where the overlinee indicate grouping of terms, and the first and second 
sec of subscripts of each element indicate the k rows and k columns, 
respectively, selected from the matrix A in order to form that particular 
kxk minor element. Expressing Che matrix A as 


A “ [a. 


1 • 


.-1 *i,'*i,+l '* \-l ’ 


V‘ -^1 "l^-" "k~* *k "k' 

the arbitrary ich column of A^ has been formed by selecting the k m- 

length columns a. , a ,..,a from A and forming all possible lock 
T. ^2 ^ 

minors from these k columns by choosing k rows at a time according to 
some predetermined convention. Thus Che arbitrary ich column of A^ is 

of length (°) , Che number of combinations of k items chosen at a time 
from m. 

Based on Laplace's Theorem [9], a possible construction for T would 
be CO make the ith row of T consist of cofactors corresponding to Che 
lock minor elements in the arbitrary ith column of A^. Thus the ele- 
ments of Che ith row of T are formed from the (m-k) complementary 
columns a^,,.,a^_j^, 

selected from these columns to form each element of the ith row T are 
complementary to the k rows selected from a. ,..,a to form the cor- 

h \ 

responding element in the ith column of A. Hence the order in which 
the (m-k) rows of A are selected in forming the elements of T may be 
said to follow a complementary convention to that used in selecting k 


T - 


'Scfl. .ffl^k+l. .m 
• 

... 

1+1.. 




o. * 

.m,l. .m-k 

• • « 




the "tilde" symbol denoting that Che minors are appropriately signed. 

Recall that we are interested in the product T B^. The matrix 

is formed from the columns of the matrix B in the same manner as A^ was 

from the columns of A. B consists of p columns of length o, as 

B • [b, . . b . , b. b. b. , b. b. , b ] , 

1 Ji-i Ji-n jfc-i Jk " 
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and Is givan as 

“b. 


b: 




1< .k,j 4 j 


4. B 


1. .k,p-kfl4 4? 


*lc tn*lc^X #• ni"*lc^ L « • s y X « • p 

Hara Che arbitrary jch coluam of B^ is formed by picking the columns 
b. , 44 ,b. from B and forming kxk minors using the same convention for 


'j 


i 


Che order In which k rows are chosen from m as before. Now consider a 
hybrid matrix 


X 


IJ 






• a ] 


and assume that we want to compute its determinant by invoking Laplace's 
Theorem. We may well select the k columns b. ,..b from X. and 

^1 ^k ^ 

from kxk minors to be multiplied by the corresponding (m-k) x (m-k) co- 

’^l^-l’^i^^+l’'*’* 


factors from the remaining columns aj^,..,a^ -1**1 +1’ 


m 


of X ^^4 But Che former kxk minors are identical to the elements in the 

arbitrary jth column of B^, and the latter (m-k) x (m-k) cofactors 

are identical to the elements in the arbitrary ith row of T. The order 
of selection of k rows to form the former elements and of the (m-k) 
rows CO form the latter elements are precisely complementary by convention, 
and thus Che product of the ith row of T with the Jch column of 

chat is, the ijch element of T B,'', is Just dec X.., that is, 

* ij 




"jl ‘ 




'jk ■ V' ■■■■ *"■ 


But the above is precisely the arbitrary numerator of the kch kind 


N 


^1 

h" ^k 

h” h 


so chat we have r* 


-1 ^ ®k 

'' ''k dec A 


^l-.^k 
“l..\ 


N 


•^1 

N 


1”\ 


“l..“k 


N 


,^m— k+1 . . ^m 


"l..“k 


^1 ^k 
u u 

. Jr- Jk 

N 

4 Jl” ^k 

' y y 

m-k+l. . m 
N 

■’k 


/l..\ 
"p-k+l..“p 
•^i ..^i 

h 

p-k+X , • p 


N 


.^m-k+l . 4 ^m 

Vk+l..“p 


/dec A. 


Hence the kch exterior map, represented in matrix form, has the 
numerators of the kth kind as Its elements. 

Conclusion 


Whereas numerators of different kinds are used in closed-loop feed- 
back design methods in industry via well-tested software packages, not 
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ouch is avsllsblc In ch« litarsturs about thalr natura othar than chat thay 
ara datarmlnants involving columns of A and B; wa hava saan in tha a- 
bova sactlon, howavar, chat thay ara, in fact, the alamants of tha appro- 
priate induced exterior map. In earlier sections wa saw the connection be- 
tween feedback loop closures and numerators of different kinds; thus, as 
the feedback loops are closed one by one, the different Induced exterior 
maps come into play. Current multivariable feedback design methods attempt 
to shape multioutput responses based only on the last exterior map and, con- 
sistent with this, present definitions of multivariable zeros involving only 
minors of the largest order In the transfer matrix. However, because the 
intermediate exterior maps are also involved when we close loops, it would 
seem that 'we would be making better use of the available information in 
shaping multi-output responses if we defined intermediate zeros from all 
the maps when we try to set up a valid multivariable counterpart of the 
familiar SISO zeros. 
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